PATH INTEGRAL MONTE CARLO APPROACH TO THE U(l) LATTICE 
GAUGE THEORY IN (2+1) DIMENSIONS 
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Path Integral Monte Carlo simulations have been performed for U(l) lattice gauge theory in (2+1) 
dimensions on anisotropic lattices. We extract the static quark potential, the string tension and the 
low-lying "glueball" spectrum. The Euclidean string tension and mass gap decrease exponentially 
at weak coupling in excellent agreement with the predictions of Polyakov and Gopfert and Mack, 
but their magnitudes are five times bigger than predicted. Extrapolations are made to the extreme 
anisotropic or Hamiltonian limit, and comparisons are made with previous estimates obtained in 
the Hamiltonian formulation. 
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I. INTRODUCTION 

Classical Monte Carlo simulations Qof the path inte- 
gral in Euclidean lattice gauge theory [2| have been very 
successful, and this is currently the preferred method 
for ah initio calculations in quantum chromodynamics 
(QCD) in the low energy regime. Monte Carlo ap- 
proaches to the Hamiltonian version of QCD propounded 
by Kogut and SusskindQ have been less successful, how- 
ever, and lag at least ten years behind the Euclidean cal- 
culations. Our aim in this paper is to see whether useful 
results can be obtained for the Hamiltonian version by 
using the standard Euclidean Monte Carlo methods for 
anisotropic lattices and extrapolating to the Hamil- 
tonian limit in which the time variable becomes contin- 
uous, i.e. the lattice spacing in the time direction goes 
to zero. The Hamiltonian version of lattice gauge the- 
ory is less popular than the Euclidean version, but is still 
worthy of study. It can provide a valuable check of the 
universality of the Euclidean results 0, and it allows the 
application of many techniques imported from quantum 
many-body theory and condensed matter physics, such 
as strong coupling expansions the t-expansion0, the 
coupled-cluster method^], the plaquette expansion 0, 
loop representation method |lfj| . and more recently the 
density matrix renormalization group 11] (DMRG). None 
of these techniques has proved as useful as Monte Carlo 
in (3+1) dimensions; but in lower dimensions they are 
more competitive. 

A number of Quantum Monte Carlo methods have 
been applied to Hamiltonian lattice gauge theory in the 
past, with somewhat mixed results. A "Projector Monte 
Carlo" approach ^^El using a strong coupling represen- 
tation for the gauge fields runs into difficulties for non- 
Abelian models, in that it requires Clebsch-Gordan coef- 
ficients for SU(3) which are not even known at high or- 
ders; and furthermore a version of the "minus sign prob- 
lem" rears its head 141. A Greens Function Monte Carlo 
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approach was pioneered by Chin et al 15] and Heys and 
Stump 16], which uses a weak coupling representation for 
the gauge fields. This approach can be used successfully 
for non-Abelian theories, and obtains estimates of pa- 
rameters such as the string tension and glueball masses 
from the correlation functions in a similar fashion to Eu- 
clidean techniques. Unfortunately the approach requires 
the use of a "trial wave function" to guide random walk- 
ers in the ensemble towards the preferred regions of con- 
figuration space^. This introduces a variational ele- 
ment into the procedure, in that the results may exhibit 
a systematic dependence on the trial wave function. We 
have previously explored [ill Il9l. I20I] a "forward- walking" 
technique^! [2^ for measuring the expectation values 
and the correlation functions, which should minimize this 
dependence; but calculations for the SU(3) Yang-Mills 
theory in (3+1) dimensions still showed an unaccept- 
ably str ong sensitivity to the parameters of the trial wave 
function [201] . For this reason, we are forced to look yet 
again for an alternative approach. 

As mentioned above our aim in this paper is to use 
standard Euclidean path integral Monte Carlo techniques 
for anisotropic lattices, and see whether useful results 
can be obtained in the Hamiltonian limit. Morningstar 
and Peardon 4] showed some time ago that the use of 
anisotropic lattices can be advantageous in any case, par- 
ticularly for the measurement of glueball masses. We use 
a number of their techniques in what follows. 

As a first trial of this approach, we treat the U(l) gauge 
model in (2+l)D, which is one of the simplest models 
with dynamical gauge degrees of freedom, and has also 
been studied extensively by other means (see Section II) . 
Path integral Monte Carlo methods were applied to this 
model a long time ago by Hey and collaborators p3l l24| . 
but the techniques used at that time were not very so- 
phisticated, and the results were rather qualitative. Very 
little has been done since then using this approach on this 
particular model, apart from a calculation by Irback and 
Peterson 

The rest of this paper is organised as follows. In sec- 
tion II we discuss the U(l) model in (2+1) dimensions in 
its lattice formulation, and outline some of the work done 
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on it previously. The details of the simulations, including 
the generation of the gauge-field configurations, the con- 
struction of the Wilson loop operators and glueball oper- 
ators, and the extraction of the potential and string ten- 
sion estimates are described in section III. In section IV 
we present our main results for the mean plaquette, static 
quark potential, string tension and glueball masses. The 
static quark potential has not previously been exhibited 
for this model, as far as we are aware. Finally we make 
an extrapolation to the Hamiltonian limit, and compar- 
isons are made with estimates obtained by other means 
in that limit. We find that indeed the PIMC method can 
give better results than other weak-coupling Monte Carlo 
methods, even in the Hamiltonian limit. Our conclusions 
are summarized in Sec. VI. 



II. THE U(l) MODEL 

Consider the isotropic Abelian U(l) lattice gauge the- 
ory in three dimensions. The theory is defined by the 
action 



S = (3 R eTr ^V 



(1) 



where 



P^(r) = [l - i?eTr{U M (r)U„(r + /i)U} t (r + fr)U£(r)}] 

(2) 

is the plaquette variable given by the product of the link 
variables taken around an elementary plaquette. The link 
variable U p (r) is defined by 



U M (r) = exp[ieaA At (r)] = exp[i0 M (r)] 



(3) 



where in the compact form of the model, 0^{r){= 
eaAfj,(r)) <E [0, 2n] represents the gauge field on the di- 
rected link r — > r + ft. The parameter (3 is related to the 
bare gauge coupling by 



= 



9 Z 



(4) 



where g 2 — ae 2 , in (2+1) dimensions. 

The lattice U(l) model in (2+1) dimensions has been 
studied by many authors, and possesses some important 
similarities with QCD (for a more extensive review, see 
for example ref. [lj])- If one takes the "naive" contin- 
uum limit at a fixed energy scale, one regains the simple 
continuum theory of non-interacting photons 26]; but if 
one renormalizes or rescales in the standard way so as to 
maintain the mass gap constant, then one obtains a con- 
fining theory of free massive bosons. Polvakovp^j showed 
that a linear potential appears between static charges 
due to instantons in the lattice theory; and Gopfert and 
Mack 28] proved that in the continuum limit the theory 
converges to a scalar free field theory of massive bosons. 



They found that in that limit the mass gap behaves as 



8n2 i n2 rim 
am D = \ —z- exp( ^v{0)) 

V 9 9 
while the string tension is bounded by 



9' 



(5) 



(6) 



where v(0) is the Coulomb potential at zero distance, and 
has a value in lattice units 



v(0) = 0.2527 



(7) 



for the isotropic case. They argue that © represents the 
true asymptotic behaviour of the string tension, where 
the constant c is equal to 8 in classical approximation. 
The theory has a non-vanishing string tension for ar- 
bitrarily large /?, similar to the behaviour expected for 
non-Abelian lattice gauge theories in four dimensions. 
For an anisotropic lattice, the gauge action becomes^] 



(8) 



where and Pa are the spatial and temporal plaquette 
variables respectively. In the classical limit 
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(9) 
(10) 



where At = at/a s is the anisotropy parameter, a s is 
the lattice spacing in the space direction, and a t is the 
temporal spacing. The above action can be written as 



S = 



r i<j 



cos Oij (r] 



—y 

At A' 



1 - cos On (r) 



(11) 



In the limit At — > 0, the time variable becomes contin- 
uous, and we obtain the Hamiltonian limit of the model 
(modulo a Wick rotation back to Minkowski space). 

The behaviour of the mass gap in the anisotropic case 
will be similar to equation JSJ Generalizing discussions 
by Banks et aZp^] and Ben-Menahempjoj . we find that 
the exponential factor takes exactly the same form in the 
anisotropic case. The only difference is that the lattice 
Coulomb potential at zero spacing for general At is 



v(0) 



d 3 k 



At 



— 7r (2^) 3 4[sin 2 (fc /2) + AT 2 (sin 2 (fci/2) + sin 2 (fc 2 /2))] 

0.2527 (At = 1) 
0.3214 (At = 0) 



(12) 
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But this result neglects the effects of monopoles with 
charges other than 0,±1 in the monopole gas, which is 
justified in the Euclidean case, but not in the Hamilto- 
nian limit [2^. [3(j • 

The Hamiltonian version of the model has been stud- 
ied by many methods: some recent studies include se- 
ries expansions |3ll. finite-lattice techniaues|33j. the t- 
expansion 0, l35l ] . and coupled-cluster techniques j^E 
I37L l38| and the plaquette expansions M as well as Quan- 
tum Monte Carlo methods [H [H I19TT39I llof . Quite ac- 
curate estimates have been obtained for the string tension 
and the mass gaps, which can be used as comparison for 
our present results. The finite-size scaling properties of 
the model can be predicted using an effective Lagrangjan 
approach combined with a weak-coupling expansion |4l| , 
and the predictions agree very well with finite- lattice data 

0. 



ably chosen intervals. We replaced approximately 30% 
of the ordinary Metropolis updates with Fourier updates 
for anisotropy At > 0.444 and 50% for highly anisotropic 
cases (At < 0.444). These moves satisfy the require- 
ments of detailed balance and ergodicity for the algo- 
rithm. 

A single sweep involves attempting N changes to ran- 
domly chosen links of the lattice, where N(= 3N t N^) is 
the total number of links on the lattice. The first sev- 
eral thousand sweeps are discarded to allow the system 
to relax to equilibrium. Figure ^ shows measurements 
of the mean plaquette for (3 = 2.0 and At = 1.0, and 
it can be seen that equilibrium is reached after about 
50,000 sweeps, with the measurements fluctuating about 
the equilibrum value thereafter. For highly anisotropic 
cases the system was much slower to equilibrate, despite 
the Fourier acceleration, and in the worst case the equi- 
libration time was of order 100,000 sweeps. 



III. METHODS 
A. Path Integral Monte Carlo algorithm 

We perform standard path integral Monte Carlo simu- 
lations on a finite lattice of size N% x N T , where N s is the 
number of lattice sites in the space direction and N T in 
the temporal direction, with spacing ratio At = at/a s . 
By varying At it is possible to change a*, while keeping 
the spacing in the spatial direction fixed. The simula- 
tions were performed on lattices with N s = 16 sites in 
each of the two spatial directions and Nt = 16 — 64 in 
the temporal direction for a range of couplings j3 = 1 — 3. 

The ensembles of field configurations were generated 
by using a Metropolis algorithm. Starting from an arbi- 
trary initial configuration of link angles, we successively 
update link angles 0^(r, t) at positions (r, t) which are 
chosen randomly each time. We propose a change A6> to 
this link angle, which is randomly drawn from a uniform 
distribution on [—A, A], where A is adjusted for each set 
of parameters to give an acceptable "hit rate" around 
70-80%. The change is accepted or rejected according to 
the standard Metropolis procedure. 

For high anisotropy (At << 1), any change in a time- 
like plaquette will produce a large change in the ac- 
tion, whereas changes to the space-like plaquettes will 
cause a much smaller change in the action. This makes 
the system very "stiff" against variations in the time- 
like plaquettes, and therefore very slow to equilibrate, 
with long autocorrelation times. To alleviate this prob- 
lem, we used a Fourier update procedure |42l |43|. Here 
proposed changes are made to space-like links which are 
designed to alter space-like plaquette values much more 
than the time-like plaquette values. At randomly chosen 
intervals and random locations, we propose a non-local 
change A0(r, r) = X sin fc(T — t ) on a "ladder" of space- 
like links extending half a wavelength (A = 27r/fc) in the 
time direction, where both k and X are randomly cho- 
sen at each update from uniform distributions in suit- 
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FIG. 1: Plot of the mean plaquette value against the number 
of configurations. 

After discarding the initial sweeps, the configurations 
were stored every 250 sweeps thereafter for later analysis. 
Ensembles of about 1,000 configurations were stored to 
measure the static quark potential and glueball masses 
at each (3 for At > 0.4, and 1,400 configurations for 
At < 0.333. Measurements made on these stored config- 
urations were grouped into 5 blocks, and then the mean 
and standard deviation of the final quantities were es- 
timated by averaging over the 'block averages', treated 
as independent measurements. Each block average thus 
comprised 50,000 - 70,000 sweeps. 



B. Interquark Potential 

The static quark-antiquark potential, V(r) for various 
spatial separations r is extracted from the expectation 
values of the Wilson loops. The timelike Wilson loops 
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are expected to behave as: 

W(r, t) ~ Z(r)exp [— TV(r)]+(excited state contributions) 

(13) 

We have averaged only over loops (xo,ro) — > 
(x + r, t ) -> (x + r, t + r) -> (x , r + r) ^ (x , r ) 
which follow either two sides of a rectangle between x 
and x + r or a single-step 'staircase' route, to estimate 
W(r, r). To suppress the excited state contributions, a 
simple APE smearing technique |44l |45| was used on the 
space-like variables. In this technique an iterative smear- 
ing procedure is used to construct Wilson loop (and glue- 
ball) operators with a very high degree of overlap with 
the lowest-lying state. In our single-link smoothing pro- 
cedure, we replace every space-like link variable by 



Ui 



P 



aUi 



(1-a) 



U x 



(14) 



where the sum over "s" refers to the "staples" , or 3-link 
paths bracketing the given link on either side in the spa- 
tial plane, and P denotes a projection onto the group 
U(l), achieved by renormalizing the magnitude to unity. 
We used a smearing parameter a = 0.7 and up to ten 
iterations of the smearing process. 

To further reduce the statistical errors, the timelike 
Wilson loops were constructed from "thermally aver- 
aged" temporal links ptjj]. That is, the temporal links 
Ut in each Wilson loop were replaced by their thermal 



averages 



Ut 



dUUexp(-S[U])/ / dUexp(-S[U}) 



(15) 



where the integration is done over the one link only, and 
depends on the neighbouring links. For the U(l) model, 
the result can easily be computed in terms of Bessel func- 
tions involving the 'staples' adjacent to the link in ques- 
tion. This was done for all temporal links except those 
adjacent to the spatial legs of the loop, which are not 
'independent'^^]. The procedure has a dramatic effect 
in reducing the statistical noise, by up to an order of 
magnitude 3 , worth a factor of 100 in Monte Carlo run- 
time. 

The Wilson loop values W(t,t) are expected to de- 
crease exponentially with Euclidean time r. A typi- 
cal plot of the logrithmic ratios of successive loop val- 
ues is shown in Figure [3 for (3 — 2.0, At = 1.0 and 
R — ||r|j =4. It can be seen that with the heavy smear- 
ing we have used, a flat 'plateau' is attained virtually 
straight away. The Wilson loops are therefore fitted with 
the simple form 

W(r, t) = ae- ry(r) (16) 
to determine the 'effective' potential V(r). 



C. Glueball masses 

Estimates for the glueball masses were obtained from 
the time-like correlations between spatial Wilson loop op- 




FIG. 2: Logarithmic ratio of Wilson loops as a function of 
r for fixed R = 4 at j3 = 2.0 and At = 1.0. The dashed 
horizontal line indicates the plateau value. 



erators $j(r)), 



C(r) = (Mr) j Mr)) 



(17) 



in a standard fashion. As the temporal separation be- 
comes large, the above correlator tends to be dominated 
by the lowest energy state carrying the quantum num- 
bers of If these quantum numbers coincide with those 
of the vacuum state, one then looks at the next higher 
energy state. So before taking the large Euclidean time 
limit, one subtracts the vacuum contribution from the 
correlator. Thus 



Mr) = Mr)- < 0\Mr)\0 > 



(18) 



is a gauge invariant, translationally invariant, vacuum- 
subtracted operator capable of creating a glueball out of 
the vacuum. As a function of the temporal separation r, 
and with periodic boundary conditions, the correlation 
function is expected to behave as 



C(t) 



<ft( T )f i (0)> 

ci(exp(-m l (r)) + exp(-m l (T - r))) 
+(excited state contributions) (19) 



where m, is the mass of the lowest glueball state in that 
sector, and T = N T a t is the extent of the periodic lattice 
in the time direction. We project out states with momen- 
tum k = and spin J = by summing over all lattice 
translations and rotations of the operators involved in 
<I>i. In the present case, we study only the lowest-lying 
'antisymmetric' (PC — - -) and 'symmetric' (PC = ++) 
glueball states, corresponding to operators <&i which are 
the sine and cosine, respectively, of the sum of the link 
angles around the Wilson loop in question. 

The statistical fluctuations in C(r) are given by 01 



C(0) 



(20) 
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Thus, the signal-to-noise ratio collapses as C(r) falls ex- 
ponentially fast with t. Hence, it becomes important to 
use a glueball operator for which the overlap with the 
glueball state of interest is strong for small lattice spac- 
ing, and such that C(r) attains its asymptotic form as 
quickly as possible. For such an operator, the signal-to- 
noise ratio is also optimal |47|. Such operators can be 
constructed by exploiting link smearing and variational 
techniques 0, H]| . 

In the strong coupling limit = 0,thc plaquette oper- 
ator U p will create a symmetric glueball state from the 
vacuum. For large f3, however, the glueball wave func- 
tions are expected to spread out and become more dif- 
fuse. To obtain a good overlap with the ground state in 
each sector at weak coupling, we need large, smooth op- 
erators (f>i on the lattice scale. An optimized operator is 
found by a variational tec hniq ue, following Morningstar 
and PeardonQ and Teper |48|. First, we calculate the 
correlation functions for square n x n Wilson loops with 
m smearings, and determine the values of m and n for 
which the ratio (C(l)/C(0)) nm is a maximum. In the 
second pass, an optimized glueball operator was found 
as a linear combination of the basic operators 4>i , 



Finally the optimized correlation function was fitted with 
the simple form 



Ci(r) = ci coshTOi(— - t) 



(26) 



to determine the glueball mass estimates. Figure [3] shows 
an example of the correlation function and fit for the 
antisymmetric state at (3 — 2 and At = 1. It can be seen 
that the form l|26|l fits the data very well. 



C(X) 



$(t) = y^^q0ia(r) 



(21) 




where the index a runs over the rectangular Wilson loops 
with dimensions l x = [n — 1,71 + 1], l y = [n — l,n+ 1] and 
smearings n s = [m— 1, m+1] with n and m as determined 
in the first pass, making 27 operators in all. The 27 x 27 
correlation matrix was computed 



Ciafi(r) = /]< 0\<t>ia(T + T )(j>ifi(To)\ > 



(22) 



FIG. 3: Antisymmetric glueball correlation function C(t) for 
the channel against r at At — 1.0 and f3 = 2.0. The solid 
curve is a fit to the simulation results using Eq. I|26|l . 



IV. SIMULATION RESULTS AT FINITE 
LATTICE SIZE 



where 4>i a (r) is a vacuum-subtracted operator 

<Mt) = <j> ia {r)- < 0|& o (t)|0 > 

The coefficients V{ a were then determined by minimizing 
the effective mass at t = 1 



772(1) 
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■In 



Eo/3 ViaVipCiapil) 



(23) 



Let Vj denote a column vector whose elements are the 
optimal values of the coefficients u, a , then the column 
vector Vj formed from the Vi a is the solution of an eigen- 
value equation 



C(l)v, 



-o t m(l) 



C(0)v 4 



(24) 



The eigenvector Vo corresponding to the largest eigen- 
value e~ atm then yields the coefficients for the operator 
$i(r) which best overlaps the lowest- lying state. 

A third pass was made to estimate the optimized cor- 
relation function 



Ci(r) = V < 0|$ 4 (r + t )$ 4 (t )|0 > 



TO 



(25) 



Simulations were carried out for lattices of Nf x N t 
sites, with N s = 16 and Nt ranging from 16 to 48 
sites, with periodic boundary conditions. Each run in- 
volved 250,000 sweeps (350,000 for high anisotropy) of 
the lattice, with 50,000 sweeps (100,000 high anisotropy) 
discarded to allow for equilibrum, and configurations 
recorded every 250 sweeps thereafter. Coupling values 
from j3 = 1.0 to 3.0 were explored at anisotropies At 
ranging from 1 to 1/3. We fixed At = 16 /N t in the 
first pass, so that the lattice size remains fixed at 16a s 
in all directions. At strong couplings (small /3), we ex- 
pect the behaviour to be the same as in the bulk system, 
but at weaker couplings (large 0) the finite-size/finite- 
temperature corrections will become more important. 
We shall monitor our data for signs of these effects. The 
results for the Euclidean case At = 1 are listed in Table 

Q 



A. Mean Plaquette 

Figure 01 shows the behaviour of the mean spatial pla- 
quette < P > for different j3, at fixed At = 1 (Euclidean, 
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TABLE I: Monte Carlo estimates for the mean plaquette < 
P >, string tension K, symmetric and antisymmetric glueball 
masses M ++ , M — , and the mass ratio Rm in the Euclidean 



case At 


= 1.0. 










P 


< P > 


K 


M 0+ + 


M - 


Rm 


1.0 


0.475 


0.674 


2.69(1) 


2.6(1) 


1.00(6) 


1.35 


0.629 


0.343(6) 


2.14(5) 


1.65(5) 


1.29(5) 


1.41 


0.656 


0.286(5) 


2.08(5) 


1.58(3) 


1.31(7) 


1.55 


0.704 


0.200(1) 


1.79(5) 


1.26(3) 


1.41(7) 


1.70 


0.748 


0.122 


1.41(3) 


0.88(1) 


1.60(8) 


1.90 


0.790 


0.082 


1.14(3) 


0.54(1) 


2.0(1) 


2.0 


0.806 


0.050 


0.79(1) 


0.44(1) 


1.78(9) 


2.25 


0.834 


0.022 


0.50(2) 


0.236(9) 


2.1(1) 


2.5 


0.854 


0.012 


0.34(2) 


0.165(9) 


2.0(1) 


2.75 


0.869 


0.009 








3.0 


0.881 


0.010 









isotropic case). A strong coupling perturbation series ex- 
pansion has been obtained for this quantity by Bhanot 
and Creutz^^ to order (3, and a weak coupling series to 
order 1//3 5 by Horsley and Wolff 50]. These series are 
represented by solid and dashed lines on the graph, re- 
spectively. It can be seen that the data follow the strong- 
coupling expansion for (3 < 1.5, and match the weak- 
coupling expansion quite closely beyond (3 ~ 2. The 
variation of < P > with coupling is extremely smooth, 
with no sign of any phase transition, as we should ex- 
pect. The cross-over from strong to weak coupling seems 
to take place quite rapidly in the region /3 « 1.8 — 2.0. 
Horsley and Wolff |H3 investigated the effect of a finite- 
size lattice in their weak-coupling expansion calculations. 
They found that such effects enter at order 1 / (3 2 as a cor- 
rection of order 1/L D , where L is the lattice size and D 
is the number of dimensions. Thus for the compact U(l) 
model in 3-dimensions and L > 10, the finite size effects 
should be essentially negligible. 

Figure [5] shows a plot of our estimates of < P > as a 
function of anisotropy jHJ 1 At 2 for the case = It 
can be seen that < P > remains almost constant. We 
would like to make contact with previous Hamiltonian 
studies by showing that the mean plaquette value ap- 
proaches to previously known values in the Hamiltonian 
limit At — > 0. The extrapolation was performed using a 
simple cubic fit in powers of At 2 . In this limit our results 
agree very well with the Hamiltonian estimate obtained 
by Hamer et aljl9j. Our general estimates in the Hamil- 
tonian limit At = are listed in Table ITU 



1 The physical anisotropy, as measured spatial versus tempo- 
ral correlation lengths for instance, will differ from the 'bare' 
anisotropy At. Since we are only interested in the Hamiltonian 
limit At — > 0, however, this is of no concern for our present 
purposes 




P 



FIG. 4: The mean plaquette as a function of f3 at At = 1. 
The circles are our Monte Carlo estimates. The solid curve 
represents the 0(/3) strong-coupling expansion |49|l and the 
dashed curve represents the 0(l//3 4 ) weak-coupling expansion 

M. 



<p> 




FIG. 5: The mean plaquette as a function of At 2 at (3 = \/2. 
The solid curve is a cubic fit to the data, extrapolated to the 
Hamiltonian limit. The triangle shows the Hamiltonian series 
estimate |l9|l in that limit. 



FigureElgraphs the resulting estimates of < P > in the 
Hamiltonian limit At = as a function of coupling (3. 
The weak-coupling an( A strong-coupling 3l| series pre- 
dictions are shown as dashed and solid lines on the graph 
respectively, while some previous Greens Function Monte 
Carlo estimates are shown as triangles. Our present 
results are generally in reasonable agreement with the 
earlier ones, if perhaps a little low in places. It can be 
seen that the crossover from strong to weak coupling be- 
haviour again occurs at around (3 ~ 1.8. 
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TABLE II: Monte Carlo estimates for the mean plaquette 
< P >, string tension K, symmetric and antisymmetric glue- 
ball masses M ++, M -~, and the mass ratio Rm in the 
Hamiltonian limit At — 0. 



p 


<P> 


K 


M 0++ 


M - 


Rm 


1.0 


0.397(6) 


0.302(3) 


1.9(1) 


1.45(7) 


1.32(1) 


1.35 


0.599(2) 


0.132(8) 


1.4(1) 


0.9(1) 


1.5(2) 


1.41 


0.628(1) 
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FIG. 6: Mean plaquette estimates as a function of /3 at 
At = 0. Our present Monte Carlo estimates are shown by 
circles. The strong-coupling [3l| and weak-coupling |4l| se- 
ries predictions are shown as solid and dashed lines respec- 
tively. The GFMC estimates [Tfll ] are represented by the solid 
triangles. 



B. Static quark potential and string tension 

Figure shows a graph of the static quark potential 
V(R) as a function of radius R at j3 = 2.0 and At = 1.0. 
To extract the string tension, the curve is fitted by a form 



V(R) =a + bh\R + aR, 



(27) 



including a logarithmic Coulomb term as expected for 
classical QED in (2+1) dimensions which dominates the 
behaviour at small distances, and a linear term as pre- 
dicted by PolvakovpT] and Gopfert and Mack 28] dom- 
inating the behaviour at large distances. The linear be- 
haviour at large distances is very clear, but the data do 
not extend to very small distances, so there is no real test 
of the presumed logarithmic behaviour in this regime. 



V(R) 




FIG. 7: The static-quark potential V(R) as a function of 
the separation R. This plot involves measurements at /? = 2.0 
for At = 1.0 with 10 smearing sweeps at smearing parameter 
a = 0.7. The errors are smaller than the symbols. The dashed 
line is a fit to the form V(R) — a + bR + cln(i?). 



Figure |H1 shows the behaviour of the fitted value of 
the string tension K = era 2 as a function of j3 for the 
Euclidean, isotropic case (At = 1). The solid line on 
the graph represents the form © predicted by Gopfert 
and Mack|28j]. using a value of c = 44.0 ± 0.4. It can be 
seen that this form represents the data rather well over a 
range 1.41 < (3 < 2.5; in fact an unconstrainted fit to the 
data gives a slope of 2.49 ± 0.15, extremely close to the 
predicted value 2.494. The coefficient c, however, is much 
bigger than the value c = 8 predicted in the classical 
approximation. It would be interesting to explore how 
higher-order quantum corrections affect the prediction 
for c. 




FIG. 8: MKyTp) as a function of j3 at At = 1.0. The 
solid curve represents the predicted asymptotic form, eq. 
with our estimated value for the normalization constant, c = 
44.0 ± 0.4. The dashed line represents the finite size scaling 
behaviour|4l|. The solid triangles show the previous estimates 
of Irback and Peterson |2a| . 
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The dashed line in Figure |S1 gives some idea of the 
expected finite-size scaling corrections to the string ten- 
sion. These have not been calculated explicitly for the 
Euclidean model, as far as we are aware, but in the 
Hamiltonian version the string tension at weak coupling 
is found j^U to behave as 



where L = N s is the lattice size (here L = 16), and this 
is represented by the dashed line. This would predict 
that the string tension will be dominated by finite-size 
corrections beyond j3 ~ 2.5, and indeed the Monte Carlo 
estimates do flatten out beyond that point, although at 
a level below equation 

Figure[!|]shows the behaviour of the string tension K as 
a function of the anisotropy At 2 , for fixed coupling (3 — 
\/2. An extrapolation to the Hamiltonian limit At — ► is 
performed by a simple cubic fit. Again the extrapolation 
agrees well with earlier Hamiltonian estimates|14j. Note 
that this quantity depends rather strongly on At: there 
is a factor of three difference between the values at At = 
and At = 1. Extrapolating to At — > 0, estimates of 
the string tension in the Hamiltonian limit are obtained 
for various (3 values. 



the string tension is consistent, within errors, with the 
finite-size behaviour predicted by equation 1)28(1 , which is 
shown as a dashed line. 




1 1.2 1.4 1.6 1.8 2 2.2 2.4 2.6 2.8 



P 

FIG. 10: Graph showing estimates of the string tension in the 
Hamiltonian limit as a function of /3. The solid line is a least 
square fit to the form K = (3~ 1 ^ 2 exp(ao — ai/3), with ao = 
2.359 and ai — 3.159. The dashed line represents the finite 
size scaling behaviour |4ll. Earlier results from an exact linked 
cluster expansion[5l| and quantum Monte Carlo simulations 
flU are shown as solid triangles and squares respectively. 
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AT 2 

FIG. 9: String tension, K, as a function of At 2 at (3 = \/2. 
An extrapolation to the Hamiltonian limit is performed by a 
cubic fit and shown by the dashed line. Our MC estimates are 
shown by triangles and an earlier Hamiltonian series estimate 
[3l| is shown by a solid circle. 

Our estimates of the string tension in the Hamilto- 
nian limit are graphed in Figure I10I together with ear- 
lier results from an 'exact linked cluster expansion' |5lj 
and a quantum Monte Carlo simulation^} . It can be 
seen that our values are consistent with earlier results, 
though less accurate, and extend further into the weak- 
coupling region. The solid line in the graph represents a 
least-square fit of the weak-coupling asymptotic form ijtjjl. 
with v(0) = 0.3214 and c = 46.7 ± 0.4. This form rep- 
resents the data well for 1.35 < (3 < 2.0. Beyond (3 = 2 



C. Glueball Masses 

Figure ^2 shows results for the antisymmetric 0~~ 
glueball mass against (3 for the isotropic Euclidean case 
At = 1. The solid line on the graph is a fit to the 
data over the range 1.4 < (3 < 2.25 using the predicted 
asymptotic form, equation (0, but with an additional 
multiplying constant: 

M = am D = ci V / 87r 2 /3exp(-7r 2 /3t)(0)) (29) 

where c\ = 5.23 ± 0.11 when adjusted to fit the data. 
Thus the slope, 2.48 ±0.09, of the data matches the pre- 
dicted asymptotic form very nicely, but the coefficient is 
too large by a factor of 5.2. It would again be interest- 
ing to explore whether this discrepancy could be due to 
quantum corrections. 

The expected finite-size scaling behaviour of the mass 
gap near the continuum critical point in this model is 
not known; but Weigel and JankejH^ have performed 
a Monte Carlo simulation for an 0(2) spin model in 
three dimensions which should lie in the same univer- 
sality class, obtaining 

M - 1.3218/L (30) 

for the magnetic gap. The dashed line in FigurelTTlshows 
this prediction for L = 16. It can be seen that the Eu- 
clidean mass gap should not be affected by finite-size cor- 
rections until (3 > 2.8. 
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FIG. 11: The scaling behaviour of the antisymmetric mass 
gap against /3 at At — 1.0. The solid line is a fit of the 
form eq. 11291 . The errors are smaller than the symbols. The 
dashed line shows the finite size scaling behaviour [H^l . 



To check the consistency of our method, we plot the 
dimensionless ratio of the antisymmetric mass gap over 
the square root of the string tension against (3 together 
with the results of Teper[48j in Figure ^] The agree- 
ment is excellent. The solid line gives the ratio of the 
fits in Figures 00 and ^JJ and shows how this ratio van- 
ishes exponentially in the weak-coupling limit, whereas in 
four-dimensional confining theories it goes to a constant. 




FIG. 12: The dimensionless ratio Mj\fK as a function (3. 
Our estimates are shown by circles and solid triangles show 
the earlier results of Teper |48|. The solid curve represents 
the predicted weak-coupling behaviour. 

Figure ^5] shows the behaviour of the glueball masses 
as functions of At 2 for [3 — \/2. The extrapolation to 
the Hamiltonian limit is performed by a simple cubic fit 
in powers of At 2 . In this limit we reproduce the earlier 
estimates of Hamer et al|3lj for the 0~~ and ++ states. 

Estimates of the antisymmetric mass gap in the Hamil- 
tonian limit At = are graphed against f3 in Fig- 



ure El Also shown are results from previous strong- 
coupling series ex trap olations [3 1 ) and quantum Monte 
Carlo calculations [U- It can be seen that our present 
results agree with previous estimates but are less accu- 
rate. The solid line is a fit to our data for 1.4 < (3 < 2.25 
of the form {2g|, with v(0) = 0.3214 and ci = 5.5 ± 0.2, 
which is similar to the coefficient found in the Euclidean 
case. The fit to the data gives a slope of 3.1 ± 0.2 and 
an intercept of 3.6 ± 0.3 of the scaling curve. We note 
that the exponential slope in previous studies is gener- 
ally somewhat less than this, as tabled in and as illus- 
trated by the black triangles in Figure ^] The dashed 
line represents the finite-size scaling behaviour, equation 
(|30|l . which we assume holds in the Hamiltonian limit 
also, for want of better information. It can be seen that 
the finite-size corrections are predicted to dominate for 
(3 > 2.2, but the date are not accurate enough at weak 
couplings to establish whether this is really the case. 




FIG. 13: Estimates of the masses of ++ and glueballs 
against At 2 . Results at (3 = \pl for the ++ and are 
labeled by circles and triangle respectively. The solid and 
dashed curves are the cubic fits to the data extrapolated to 
the Hamiltonian limit. The series estimates of Hamer et al 
[3H ] in the limit At — > 0, for symmetric and antisymmetric 
channels are shown as a star and diamond respectively. 

Finally, Figure IT51 displays the behaviour of the dimen- 
sionless mass ratio, 



Rm = Af(0++)/M(0~) 



(31) 



for the Euclidean case Ar = 1. As in the (3+l)D confin- 
ing theories, we may expect that quantities of this sort 
will approach their weak-coupling or continuum limits 
with corrections of 0(l/a e //), where a e ff is the effective 
lattice spacing in 'physical' units when the mass gap has 
been renormalized to a constant. Hence for our present 
purposes we define a e ff from equation J^J) as 



a eff = v / 87T 2 /3exp(-7r 2 ^(0)) 



(32) 



with v(0) = 0.2527 for the Euclidean case. The mass 
ratio is plotted against a e // in Figure [TBI At weak cou- 
pling, we expect the theory to approach a theory of free 
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FIG. 14: Hamiltonian estimates of the antisymmetric mass 
gap plotted as a function of 8. The solid curve is the fit to the 
data for 1.4 < 3 < 2.25. The dashed line represents the finite 
size effects |5^j. The previous results from series expansion 
(I^H and quantum Monte Carlo calculations |l4jl are shown as 
solid triangles and open squares respectively. 



bosons 28] so that the symmetric state will be composed 
of two 0~~ bosons and the mass ratio should approach 
two. Our results show that as a e // goes to zero, the mass 
ratio rises to around the expected value of 2.0. A linear 
fit to the data from 0.08 < a e // < 0.32 gives an inter- 
cept Rm — 1-95 ± 0.05. However, we note that the last 
two of our estimates, together with two from Teper[4^|. 
lie considerably above Rm = 2. In the bulk system, of 
course, the ratio cannot rise above 2, because it is al- 
ways possible to construct a ++ state out of two 0~~ 
mesons. These points correspond to couplings (3 > 2, 
and we conjecture that they may be affected by finite- 
size corrections: a simulation on a larger lattice would 
be necessary to check on this point. Our last two points 
have not been included in the fit. 

Figure ^| shows a similar graph for the Hamiltonian 
limit, At = 0. Within errors, our present results are 
consistent with earlier series [31) and quantum Monte 
Carlo estimates, but are much less accurate, and tend 
to lie consistently on the high side. A linear fit to the ear- 
lier data from 0.02 < a eff < 0.12 gives R M = 2.14±0.01. 



FIG. 15: A graph showing estimates of the mass ratio Rm 
as a function of the effective spacing, a e //, at At = 1.0. Our 
present estimates are shown by the circles. The dashed line 
is a linear fit to the data over the range 0.08 < a e ff < 0.31. 
The solid triangles show the previous estimates of Teper |48|l . 




FIG. 16: Mass ratio in the Hamiltonian limit as a function 
of the effective spacing, a e //, at At = 0. Our MC estimates 
are shown b y th e circles. Earlier series [3l| and quantum 
Monte Carlo |2Q] results are shown by solid triangles and open 
squares respectively. The dashed line is a linear fit to the 
earlier data from series expansions [3l| over the range 0.02 < 
a e ff < 0.12. 



V. CONCLUSIONS 

In this paper, we have applied standard Euclidean 
path integral Monte Carlo methods to the U(l) model 
in (2+1) dimensions on an anisotropic lattice, and taken 
the anisotropic limit At — * to obtain the Hamiltonian 
limit of the model. 

We have obtained the first clear picture of the static 
quark potential in this model, showing very clear evi- 
dence of the linear confining behaviour at large distances 
predicted by Polvakov|27|. There is also a turnover at 
short distances consistent with a logarithmic Coulomb 



behaviour in that regime. 

In the isotropic or Euclidean case At = 1, the string 
tension and mass gap display an exponential decrease 
at weak couplings which is in excellent agreement with 
the behaviour predicted by Polvakov[27j and Gopfert and 
Mack|2S||. Both quantities, however, are 5-6 times larger 
in magnitude than the theory predicts. It would be in- 
teresting to calculate whether higher-order quantum cor- 
rections can account for this discrepancy. 

The dimensionless ratio Mj\fK scales exponentially 
to zero in the weak-coupling or continuum limit, as pre- 
dicted by the theory. The mass ratio of the two lowest 
glueball states scales against the effective lattice spacing 
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towards a value close to 2.0, as expected for a theory of 
free scalar bosons, apart from some anomalous results at 
large (3 which we have ascribed to finite-lattice effects. 

In the anisotropic or Hamiltonian limit At — * 0, our 
results are less accurate, because of the extrapolation 
needed to reach this limit. Nevertheless, the results are 
generally in good agreement with those obtained by other 
methods. Once again, an exponential behaviour of the 
string tension and glueball masses can be demonstrated 
at weak coupling. The dimcnsionless mass ratio again 
scales to a value near 2.0. Because the exponential slope 
is steeper, finite-size effects seem to be somewhat more 
important in the Hamiltonian regime than in the Eu- 
clidean one. 

Our major object in this study was to compare the Eu- 
clidean PIMC approach to quantum Monte Carlo meth- 
ods such as Green's Function Monte Carlo (GFMC) for 
obtaining estimates in the Hamiltonian limit. The PIMC 
approach suffers from the disadvantage that an extrap- 
olation is necessary to reach the limit At — * 0; while 
GFMC suffers from the major disadvantage that it relies 
on a 'trial wave function'. In the event, we have obtained 
much better results using PIMC. A clear and consistent 
picture of the string tension and glueball masses was ob- 
tained at weak coupling. Using GFMC, on the other 
hand, only qualitative estimates of the string tension 
were obtained, and the glueball mass estimates were vir- 
tually worthless^^. No doubt there are many 'tricks of 
the trade', such as smearing and variational techniques, 



which could be used to improve the GFMC results; but 
we found previouslv|20j that there is an unacceptably 
strong dependence on the trial wave function using that 
technique, especially for large lattice size. The PIMC 
technique seems to offer a much more robust and un- 
biased approach to Hamiltonian lattice gauge theories. 
Of course, one must generally expect the Hamiltonian 
estimates to be less accurate than the Euclidean ones 
because of the extra extrapolation involved. 

We note that the PIMC results are still less accurate 
than some older quantum Monte Carlo results of Hamer, 
Wang and Price [lj|- The latter were obtained using 
a strong-coupling representation, however, and this ap- 
proach has been found to fail for non-Abelian models 
due to the occurrence of a 'minus sign' problem, as men- 
tioned in the introduction. 
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